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This short note describes a method to tackle the (bipartite) quantum separabihty prob- 
lem. The method can be used for solving the separability problem in an experimental 
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. setting as well as in the purely mathematical setting. The idea is to invoke the following 

, characterization of entangled states: A state is entangled if and only if there exists an en- 

pv^j \ tanglement witness that detects it. The method is basically a search for an entanglement 

I witness that detects the given state. 

^ : 

00 1 Introduction 

I ■ 
I \ 

^ \ Entangled quantum states are interesting both from theoretical and practical points of view. 

^ \ Theoretically, entanglement is connected with the confounding issue of nonlocality. Practically, 
■ entangled states are useful in quantum cryptography and other quantum information processing 
p • tasks (for example, see PP and references therein). The problem of determining whether a state 
Jj . is entangled or separable is thus important and comes in two flavors - one mathematical, and 
^ ! the other experimental. 

^ \ The mathematical problem of separability (for bipartite systems) is defined as follows. Let 

R^' 'Hm,n denote the set of all Hermitian operators mapping ® to ® C^. The set of 
^ ^ ■ bipartite separable quantum states Sm,n in 'Hm,n is defined as the convex hull of the separable 
^ . pure states ® € 'Hm,n}, where \a) is a norm-1 vector in C*^ and is a norm-1 

^ ! vector in C^. The set of separable states Sm,n may be viewed as a compact, convex subset 

of M^^^ by expressing each density operator as a real linear combination of the canonical 

Hermitian generators of SU{M) and SU{N) j^. The quantum separability problem is now 
easily defined as an instance of the Weak Membership problem [21 : Given a convex set K C M", 
a point p e M", and an accuracy parameter 6 > 0, assert either that (i) p e S{K, 6) (i.e. p 
is "almost in" K) or that (ii) p ^ S{K,—S) (i.e. p is "not 5-deep within" K), where S{K,S) 
denotes the union of all 5-balls with centers belonging to K, and S{K, —6) denotes the union 
of all centers of all (5-balls contained in K (in the standard Euclidean norm). The separability 
problem has been shown to be NP-hard |3], thus any devised test for separability is likely 
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to require a number of computing resources that scales exponentially with M and A^. There 
exist efficient "one-sided" tests for separability, where the output of some polynomial-time 
computable function of the matrix for p can indicate that p is certainly entangled IHl El E] 
or certainly separable [HI El t>ut not both. The algorithm in jH] is in principle a one-sided 
test because it requires an infinite amount of computational resources to detect some entangled 
states. 

The experimental flavor of the separability problem can be defined as follows: Given many 
physical copies of a completely unknown quantum state p G Hm,n, determine whether p is 
separable. One way to solve this problem is to perform a full state tomography in order to 
construct the density matrix for p to some precision 5 > 0, and then solve the mathematical 
separability problem. If rather there is some partial knowledge of p, then there are certainly 
more options, such as testing for a violation of a specific Bell- type inequality [121 HH] ot invoking 
entanglement witnesses [TH IT3] . As well, in the case where MN < 6, the positive partial 
transpose (PPT) test [3 E] can be implemented physically ^71 ^j, though currently this 
approach is not experimentally viable. 

In section|21we describe a method to tackle the mathematical separability problem in general. 
The idea is to invoke the following characterization ^H] of entangled states: A state p is 
entangled if and only if there exists an entanglement witness ^Hj that detects it. The method 
is basically a search for an entanglement witness that detects the given state. In section [HI 
we describe how to use this method as a novel tool for solving the experimental separability 
problem. We conclude with a brief discussion, and point to some future directions of research 
in sectional 

2 Solving Quantum Separability 

In this note we use the following definition of "entanglement witness", which differs slightly 
from the definition used in the literature. 

Definition 1. An entanglement witness is any operator A G 'Hm,n for which there exists a 
state p G 'Hm,n such that 

tr{Aa) < tr{Ap) V cr G Sm,n- 

Recalling that 'Hm,n is isomorphic to M*'^^^^, the above definition implies that for entangled 
p there exists a hyperplane in M.'^'^ ^ which separates p from the set of all separable states Sm,n- 
If we define the function 

bA '■= max tr{Aa), (1) 



2 



then the set {x G 'Hm,n '■ tr{Ax) = Ba} is one such hyperplane. The function Ba is imphcitly 
at the heart of the definition of "entanglement witness" as it pins down which, if any, of the 
hyperplanes, with normal A, separate the state p from Sm,n- The hyperplane defined by A and 
Ba is tangent to Sm,n and is thus the optimal hyperplane with normal vector A that separates 
p from Sm,n- 

The function 6^ leads naturally to an algorithm for quantum separability as follows. For A 
such that tr(yl^) = 1, define the function dA,p as 



Geometrically, dA,p is the signed distance from the state p to the hyperplane defined by fe^- 
By using equations ^ and |2l and the definition of "entanglement witness" , it follows that p is 
entangled if and only if 



Thus we have reduced quantum separability to the task of finding an A such that dA,p is neg- 
ative. By observing that the ability to calculate Ba gives an oracle for the Weak Optimization 
problem. Theorem 4.4.7 from j3] shows that it is possible to solve quantum separability with 
only polynomially many evaluations of the function Ba- Thus, the "hardness" of quantum sep- 
arability is contained in the "hardness" of evaluating Ba- In practice, and for low-dimensional 
applications, well-known sophisticated techniques (such as simulated annealing or interval anal- 
ysis) for finding global extrema would likely be sufficient to calculate Ba- Unfortunately the 
algorithm derived in [2] is not implementable on fixed precision computers, and it uses many 
more evaluations of Ba than is strictly necessary. In a later paper, we hope to demonstrate a 
better algorithm which uses far fewer evaluations of Ba and thus, hopefully, is of practical use 
for low-dimensional problems. 



In this section, we briefiy show how the above approach may be used when only partial in- 
formation about the state p G 'Hm,n is available. This is of particular use in an experimental 
setting. 

The state p can be written 



dA,p ■=Ba- tr{Ap). 



(2) 



3 A : dA,p < 0. 



(3) 



3 Solving Separability with Partial Information 




(4) 
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where pij G M and the and \^ (as defined in j2|) respectively generate the special unitary- 
groups SU{M) and SU{N). In this case, the coefficients pij are simply related to the expected 
values of Af ® Af : 

p., = (Af ® Af )/4 := tr(Af ® Af p)/4. (5) 

Let A:= {Xf®Xf}i =o,i,...,A/2_i;j=o,i,...,Ar2-i- The expected values of all elements of A constitute 
complete information about p. Suppose only partial information about p has been obtained by 
an experimental procedure, that is, only the expected values of the elements of a proper subset 
T of A are known. 

It helps to think of each density operator as a real vector of its expected values. With only 
\T\ expected values known, we now effectively project all the density operators onto span(T) 
by ignoring the components of the real vectors that correspond to the unknown expected 
values of p. Now each density operator, including our unknown p, is represented by a point 
in a |T|-dimensional "expectation space". Note that the set of points in this projective space 
representing all separable density operators is still a convex set. Call this convex set St and 
denote its elements by o" G M'"^'. Similarly, let p be the |T|-dimensional real vector of known 
expected values of p. To represent the Hermitian operator A = Xlxgr in this space, we can 
use the real vector A of the coefficients ax in order to have the correspondence tr{Ap) = A ■ p, 
where "■" denotes the standard dot-product of two real vectors. 

It is clear that the functions La and dA,p in the previous section can be redefined for the 
|T|-dimensional space: 

'■ = niax A ■ a 

and 

dA,p ■■=hA- A-p. (7) 

If there exists A such that d^^p < 0, then p is entangled; otherwise, more information is needed 
to determine the separability of p. As expected values are being gathered through experimental 
observation, they may be input to a computer program that searches for d^.p < 0. Finally, 
we point out that the idea of searching for an entanglement witness in the span of operators 
whose expected values are known was discovered independently and applied, in a special case, 
to quantum cryptographic protocols in PU] . 

4 Discussion 

By looking at quantum separability as a mathematical problem in the real Euclidean space 
]^A/ Af slightly altering the definition of entanglement witness, we have show that quan- 
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turn separability can be solved in oracle-polynomial time, for a rather natural looking oracle. 
This allows us to highlight the fact that the "hard" part of quantum separability is contained 
in the function Ba- The method we describe for solving quantum separability also gives experi- 
mentalists a tool for potentially determining if an unknown state is entangled by measuring only 
a subset of the expected values which completely describe the state. This method effectively 
trades quantum resources (additional copies of p) for classical resources (a computer able to 
calculate 6^). As well as providing a practical, implementable algorithm for low-dimensional 
quantum separability, some open questions which we hope to address in the future include 
getting a tight (exponential) bound on the complexity of calculating Ba, and determining the 
average number of expected values required to detect a random unknown quantum state. 
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